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Abstract 

The evaluation of the average number SAr(t) of distinct sites visited up to 
time thy N independent random walkers all starting from the same origin on 
an Euclidean lattice is addressed. We find that, for the nontrivial time regime 
and for large N, Sjyit) ~ Sj\f{t){l — A), where S]\f{t) is the volume of a hy- 
persphere of radius {ADtlnNy/^, A = i ^^^^ In"" ^^^^q s^^hn" In iV, d 

(n) 

is the dimension of the lattice, and the coefficients Sm depend on the dimen- 
sion and time. The first three terms of these series are calculated explicitly 
and the resulting expressions are compared with other approximations and 
with simulation results for dimensions 1, 2, and 3. Some implications of these 
results on the geometry of the set of visited sites are discussed. 

PACS numbers: 05.40.Fb,05.60.Cd,66.30.Dn 



1 



Typeset using REVT^ 



I. INTRODUCTION 



Usually, the extremely successful theory of random walks is only concerned with problems 
that involve a single {N = 1) random walker. A solid reason for this is the understanding that 
the average properties of the single diffusing walker serve to describe the global properties 
of systems formed by many walkers. However, there are other interesting diffuson problems 
that involve many random walkers for which the diffusive behavior of every walker of the 
total of is relevant, i.e., diffusion process that cannot be described by averaging over 
the properties of a single walker [|l|. The problem of evaluating the time spent by the 
first j particles out of a total of N to escape from a given region is a clear example 
Another important example, which is the subject of this paper, is the problem of evaluating 
the average number S'7v(t) of distinct sites visited (or territory explored) by a set of A^ 
independently diffusing random walkers up to time t 

The case A^ = 1 has been studied in detail since it was posed by Dvoretzky and Erdos 
^ and is discussed in many general references fiowever, the multiparticle (A^ > 1) 

version of this problem has been systematically treated only after the pioneering works of 
Larralde et al. [|,^ . These authors addressed the problem of evaluating the territory covered 
by a set of A^ independent random walkers, all initially placed at the same point, that diffuse 
with steps of finite variance on Euclidean lattices. They found asymptotic expressions for 
SN{t) for A^ ^ 1, and described the existence of three time regimes. Their results can be 
summarized as follows: 



where x = N for d = 1, x = N/ In t for d = 2 and x = N/ \/i for d = 3 . The properties 
of Si{t) are well known; in particular, Si{t) ~ t^/^ for d = 1, Si{t) ~ t/ Int for (i = 2 and 
Si{t) ~ t for d = 3. In the very short-time regime (t ^ tx), or regime I, there are so 
many particles at every site that all the nearest neighbors of the already visited sites are 
reached at the next step, so that the number of distinct sites visited grows as the volume 
of an hypersphere of radius t, S^it) ~ f^- The regime III (t^ <C t), or long-time regime, 
corresponds to the final stage in which the walkers move far away from each other so that 
their trails (almost) never overlap and S]\i{t) ~ NSi{t). The crossover time from regime I 
to regime II is given by tx ~ In A^ for every lattice. This can be easily understood if we take 
into account that the number of particles on the outer visited sites for very short times will 
decrease as N/z'^, where z is the coordination number of the lattice, so that the overlapping 
regime will break approximately when N/z^ ~ 1 or, equivalently, tx ~ InA^. Regime III 
never appears in the one- dimensional case (i.e., ~ oo), but ~ for c? = 2 and t^ = A^^ 
for d = 3. These crossover times will be obtained readily from the mathematical formalism 
discussed in the present paper. The most interesting regime is regime II (tx ^ t ^ tx), or 
the intermediate regime. For this time regime we will obtain explictly the main term and the 
first two corrective terms of the asymptotic expression of S'Ar(t) for A^ ^ 1. Higher corrective 
terms could be calculated as our method allows them to be obtained in a systematic way. 
The contribution of these corrective terms cannot be ignored even for very large values of 
A^ because they decay logarithmically with A^. However, as we will see in section 0, the use 
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of two corrective terms leads to a very good agreement with simulation results for relatively 
small values of (iV > 100). 

The paper is organized as follows. The asymptotic evaluation of S'iv(t) for a (i-dimensional 
Euclidean lattice is discussed in detail in Sec. |I[ Some geometric implications of this re- 
sult are discussed in Sec. fTT[ In Sec. |V|, we compare our zeroth- (i.e., main), first- and 
second-order term approximation for S'Ar(t) with other approximations and with computer 
simulations for one-, two- and three-dimensional simple Euclidean cubic lattices. The paper 
ends with some remarks on the applicability of this method to other diffusion problems and 
different media. Some technical details are discussed in an Appendix. 



II. THE NUMBER OF DISTINCT SITES VISITED 

We consider a group of random walkers starting from an origin site r = at time 
t = 0. A survival probability, Ti^{t, r), is defined as the probability that site r has not been 
visited by the random walkers before time t. Similarly, we can define a mortality function, 
1 — rjv(t, r), as the probability that site r has been visited by at least one walker in the time 
interval (0,t). The relationship between the number of distinct sites visited, S'Ar(t), and the 
survival probability is 

SN{t)=Y.{^-^N{t,r)} . (2) 

r 

For independent random walkers, we have Tj^(t,r) = [Tt{r)]^ , where r((r) = ri(i(:,r) is the 
one-particle survival probability. Next, the discrete analysis implicit in Eq. (|^) is replaced 
by a continuous one. Thus, we write 

dvor'^^^dr , (3) 

where vq is the volume (i.e., the number of lattice sites) of the hyphersphere with unit radius. 
It has been found for Euclidean lattices that 

r,(r) ^ f,(0 = 1 - AC'^e-'^"^' + E /^nr'") , (4) 

for ^ = r/y/2dDt ^ 1. Here D is the diffusion coefficient defined through the Einstein 
relation (r^) ^ 2dDt, t oo, with (r^) being the mean-square displacement of a single 
random walker. The values of A, fi and hi for d = 1, 2 and 3 are shown in Table |. A change 



to the new variable ^ and integration by parts (taking into account that r((oo) = 1), yields 

Sj^{t) = voi2dDtf^'^JN{d; 0, oo) , (5) 

where 

Md;a,b) = I' NTr~\0^^'d^ , (6) 



In order to evaluate the asymptotic behavior of Jjv((i; 0,oo) it is convenient to make the 
decomposition 



SN{t) 



1 - rf (r) 
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JAr((i;0,Oo) = JAr(rf;0,^x) + JAr((i;^x,oo) , 



(7) 



where is a value that should satisfy the following two conditions: (a) ^x is large enough 
so that Tt{r) can be well approximated by its asymptotic approximation Tt{C,) for > ^x, 
and (b) small enough so that 



^f(ex) = l/iv^ 

with p > 1 (say p = 2). From Eq. (^) it is straightforward to see that 

e'x -IniV 



(8) 



(9) 



satisfies both conditions. On the other hand, because at most dT/dC, = 0{1), and is a 
monotonous growing function, JNid; 0, ^x) is bounded by a term that goes as A^rf^~^(^x)^x' 
or equivalently, from Eq. (§), by a term that goes mainly as N^~p. But shortly we will show 
that JN^d; ^x, oo) goes essentially as In'^^^ A^; this means that JN^d; 0, ^x) is asymptotically 
smaller than any term in the asymptotic expansion for JN{d', 0, oo) and thus we can write 



JAr(rf; 0, oo) ^ JN{d] ^x, oo) , > 1 . 



(10) 



The previous discussion is illustrated in Fig. [1| for the one-dimensional case. In this figure 
we have represented the integrand of Jn{^'i 0, oo) for increasing values of N using as survival 
probability the exact value Tt{C) = erf(^/-\/2) and the asymptotic expression of Eq. (^ 
up to first order (n = 1). Notice that the area below the solid [broken] curve is just the exact 
[asymptotic approximate] value of SN^t) / [SDty/'^ = Jn{^', 0, oo). The value of ^x as given by 
Eq. (11) with p = 2 is marked with a symbol. It is clear from the figure that, for large N, (a) 
the integrand of JAr(l;^xOo) is well represented by the asymptotic expression of r((^), and 
(b) that, as stated for the general case in Eq. (p!OD, Jn{^', 0, ^x) Jn{^] ■CxOo) ~ Jn{^] 0, oo). 

From Eq. (H), one easily finds that 



dm) 
d^ 



-2n 



(11) 



n=0 



with jo = ji = fi, j2 = hi, . . .By inserting Eq. (0) into Eq. (||) one has the following 
expansion for Jiy{d; C,^, oo): 



JN{d; ^x, oo) ^ 2iV £ JuKn-i {d - 2n + 1) 



with 



n.=0 



? X 



By means of the substitution 

we get a more convenient expression for Ki\i{a): 



(12) 



(13) 



(14) 
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(e- - 1) 



(15) 



where, from Eq. (^), Zx ~ \nN/N. The integral in (|15]) is of Laplace type but it is not 
possible to use Watson's lemma directly to get its asymptotic behavior because C,°'{d^/dz) 
has a logarithmic singularity at 2; = ||T^. The evaluation of K]sr{a) requires the inversion 
of (|14D to obtain ^{z). By using Eqs. (|) and (p^ we get 



d 



-e + In A + /ilnr' + In 1 + ^ /i„r'" = In (l - 



(16) 



n=l 



The function ^ {z) can be readily obtained from this equation to first approximation: Notice 
that, as long as 



^2 > \ \nA\ 



(17) 



the left hand side of Eq. ([T6|) can be approximated by —d^ /2, so that the first-order solution 



to Eq. (0) is ^2 



— 21n[l — exp(z)]/(i. Equation (p!6| ) can be systematically solved in 



order to get higher-order approximations (see Appendix). The result is 

00 



(18) 



n=l 



where x = -(c?/2)/ln[l -exp(-^)]. The substitution of Eq. (|18D into Eq. (|T5|) (see Eq. (^ 
in the Appendix) yields 



KAfia) 



2(a-l)/2 



J2 J2 V^rrrfci")/ { - - n- -,m-N 



n=Q m=0 



rf(a+l)/2 



a 



where 



I{n, m; N) 



dze 



-Nz, 



In 2)" In" 



Inz). 



(19) 



(20) 



The evaluation of lN{n,m; N) for N 00 has been discussed in 0,1^. For the sake of 
completeness, we give here explicitly their expressions up to the order required to find SN{t) 
to second order in 1/ In A^: 



I{n,0-N) 
I{n,2-N) 



1 

N 
1 

N 
1 

N 



wy nfra — 1) 7^ + 7r^/6 
1 + ^47 + ^-^ - ' o J +•• 



ln"A^ 
In" AT 
(In" N)hi^hiN + 



In AT 



InlnA^ 1 + 



n7 
h^ 



In^A^ 

7 

^ h^ ^ 



(21) 

(22) 
(23) 



where 7 ~ 0.577215 is the Euler constant. Using these results we get from Eqs. (^), (|T2D and 
(p!9|) the following expansion for the average number of distinct sites visited on a Euclidean 
lattice of dimension d: 



^^(t)^^jv(t)(l-A) 



(24) 
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with 



d/2 



1 oo n 

A = A{N,t) = - ^In-'^AT ^ In^lniV 



71=1 



m=0 



and where, up to second order (n = 2), 

-duj. 



Si'^ = dfi, 



.(2) 



J2) 
-Si 



,(2) 




T2 + Tj-n^"^". 



(25) 
(26) 

(27) 
(28) 

(29) 
(30) 
(31) 



Here u = ^ + InA + filn{d/2), and A, fi and hi are given in Table | for d = 1, 2 and 3. 
Notice that the time dependence of A{N,t) comes from the term u through the function 
A{t). However, this function does not depend on time for the one-dimensional case and thus 
A only depends on A^. 

Recently, Sastry and Agmon have found an approximate formula for S^it) for the 
one-dimensional case. The straightforward method used by these authors is based on the 
fact that the function rf^(r) that appears in the integrand of Eq. (|^) approaches a step 
function when N — >■ oo. In this way they found 



SN{t) 



A\/15i\J In - In Va InA^ 

1 In In A^ — In a 



'--4 



InN 



+ 



(32) 
(33) 



where a is given by a = 7rexp(— 2/7r) 
the first-order approximation of Eq. ( 



1.66. It is instructive to compare this formula with 
for the one-dimensional case: 



SN{t) ^ AVDtlnN 



1 - 



lln InA^ 



In AT 



2uj ,^f\n'\nN 



(34) 



Note that the prefactor 4(Dt In A^)^/^ of the formula of Sastry and Agmon is in agreement 
with that of Eq. (plf). In Ref. [|1^, they found it "amusing" that the value a = 1 produces 
very good agreement between the approximation of Eq. (^) and the exact numerical inte- 
gration. Our Eq. (|3^) enlightens this point: Comparing Eqs. (|33|) and Eq. (|3l) one sees that 



In a is playing the role of 2lj. But u = •y — ^Irni = 0.0048507- ■ ■ for the one-dimensional 
lattice, so that In a when a = 1 leads to a good approximation to the rigorous coefficient 
2u). The equation of Sastry and Agmon for a = 1 and our first-order approximation should 
thus be very close. This is clearly confirmed in Fig. ^. 

A question to be answered is why Eq. (^) is valid for time regime II only, i.e., why it 
is not always valid for arbitrarily large values of time. The reason is that our formulas have 
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been obtained by assuming that the condition (0) holds for those values of ^ which are inside 
the integration interval [^xjOo] of the relevant integral J]\[{d;C,x,oo) that is responsible for 
the asymptotic behavior of S^lt). This implies that for our procedure to work, it is necessary 
that > InA or, from Eq. (|), that 

\nN^\\nA\. (35) 

Thus we can estimate the time Tx for which our method breaks down by solving | In A{t^ ) | ~ 
In A^. From the expressions for A quoted in Table | one finds that Tx ~ for d = 2 and 
Tx ~ A^^ ioY d = 3. For d = 1 and large N, the condition ( p5D always holds because 
A = (2/7r)^/^ is a constant and then Tx = oo. We see that the upper times Tx beyond which 
Eq. is no longer valid coincide with the crossover times tx defined in Sec. |, so we can 
say that the expressions for S^lt) given in this paper are valid only in the time regime II. 
This means that our procedure marks its own limit of validity as that of regime II, thus 
predicting the existence of a crossover time in a natural way, i.e., as a consequence of the 
mathematical formalism. 



III. GEOMETRIC PROPERTIES OF THE EXPLORED REGION 

In this section we will give a geometric interpretation of the main result of this paper, 
namely, Eq. (p^) . The quantity S]\f{t) is by definition the volume of the region Q explored 
by N random walkers after a time t from their initial deposition on a given site of the lattice 
(if the length of the lattice bonds is taken as the unit). For very short times (regime I or 
t <S IniV) the exploration is performed in a compact way because all the neighbor sites 
of any visited site are always visited at the next time step. Therefore, the explored region 
Q is an hypersphere whose radius grows ballistically and its volume is proportional to t^. 
After the regime II is reached, the development of two qualitatively different zones in the 
explored volume is observed: (i) an hyperspherical compact core of visited sites, and (ii) a 
corona of dendritic nature characterized by filaments created by those relatively few walkers 
that are wandering in the outer regions, i.e., wandering at distances significantly larger than 
the root-mean-square displacement (r^)^/^ = y2dDt of a single walker. Figure ^ shows a 
snapshot of the set of sites visited by A^ = 1000 random walkers at time t = 900 (every 
walker makes a jump at each time unit) for dimension two. The visited sites are in white 
and the inner black and outer white circles delimit the corona. The radius i?+ of the outer 
circle is equal to the maximum displacement from the origin reached by any of the walkers 
at time t. It has been argued in ||Q that the volume of this outer circle is on average given 
by the main term of (p^, i.e., by S]\f{t) = vo^ADtlnNY^'^. From this statement we can draw 
two conclusions: First, that the average radius of this outer circle is 

{4Dt\nN)^^\ (36) 

and second, that the asymptotic corrective terms (given by A) to S]\f{t) account for the 
number of unvisited sites that are inside the corona. In other words, A is the fraction of the 
volume inside the external circumference that has not been visited by any of the A^ random 
walkers. This result can be used [|l^] to easily estimate that the thickness of the dendritic 



corona is approximately given by -R+A. 
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It is also noteworthy that A{N,t) depends on t very smoothly in the time regime II as 
this dependence is due to terms proportional to powers of\nA(t) [and A(t) does not change 
exponentially: see Table ^. For the two-dimensional case, this statement is especially valid 
because A{t) ~ 1/lnt. Therefore, the ratio (given by A) between the radial size of the 
corona of Q and the radial size of fl itself remains almost constant throughout the time 
regime II. This implies that a conveniently scaled sequence of snapshots of the set of visited 
sites should be very similar (in a statistical sense), i.e., we find that fl grows, to a large 
extent, in a self-similar way inside time regime II. This property is illustrated in Fig. |^. 
As Eq. ( |5BD shows, the appropriate scale factor must be proportional to This "almost" 
self-similar behavior disappears as the regime III is approached because the correction to the 
main term of S'at (t) becomes as large as this main term, i.e., because A{N, t) approaches the 
value 1. This transition takes place when t ^ as follows from (p6D, i.e., this value coincides 
with the threshold for regime III deduced in the previous section. From the geometric point 
of view this transition corresponds to the breaking of the self-similar growing behavior by 
the appearance of a corona of filaments as large as the compact core, which finally gives 
rise to a set of separated trails that (almost) never more overlap. For the two-dimensional 
case the transition time from regime II to regime III is so great for any significant number 
of walkers that it can not be studied by numerical simulation. 



IV. NUMERICAL RESULTS 

We carried out numerical simulations for the number of distinct sites visited by = 2™, 
with m = 0, 1, . . . , 14 in two and three dimensions. For the one- dimensional case it is not 
necessary to carry out simulations because the survival probability is exactly known on this 
lattice, Tt{r) = erf (^C,/V^^ , and therefore the integral for 3^(1) as given by Eq. (|^) can be 
computed numerically. 

In our simulations, the random walkers are placed initially at the center of an hypercubic 
box of side L. The regime II is reached almost immediately with the number of random 
walkers we have used {tx ~ 10 for A^ = 2^^). The simulations were carried out only to a 
maximum time t = 200 which is sufficient for the stabilization of regime II conditions. The 
square box side for d = 2 was taken to be L = 400 to avoid any random walker reaching 
the edge of the box before the maximum time t = 200. Memory limitations forced us to 
reduce the box side to L = 200 for the three-dimensional case. While this implies a possible 
appearance of finite-size effects, we can consider them as to be negligible because the average 
displacement of the random walkers at the maximum time is small compared with L/2. Each 
experiment was repeated 10'^ times in order to achieve reasonable statistics. 

Results are plotted in Fig. ^ for one, two and three dimensions. The dots are the 
simulation results (numerical integration results in one dimension) and the broken and solid 
lines are the prediction of Eq. (E3) to first and second order, respectively. The crosses 



are the results of Sastry and Agmon |TT[ given by Eq. (|3^) with a = 1. The dotted lines 



correspond to the result of Larralde et al. given by Eq. (|I|) using the correct amplitude of 



the main term (see |jT2|). The quantity plotted is 

-1 2/d 



N 



s 



N 



(37) 
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versus 1/ In A^. From Eq. (p^) one sees that the theoretical prediction for this quantity 
is S ~ {l/d){l — A)^/'^. The agreement between the second-order approximation and the 
simulations is found to be excellent for N ^ 100. Good agreement for lower values of N 
would be expected if higher-order terms in the series were included. The importance of the 
corrective terms is evident. For example, for the one-dimensional case, we would need to 
use values of as large as 10^^ in order to obtain the same precision with the main term 
as we get with the main and two corrective terms for values of as small as 2^. Similar 
statements can be made for the other lattices, as Fig. ^ shows. 



V. REMARKS 

In this paper we have developed a method for calculating the mean number of distinct 
sites visited by independent random walkers on Euclidean lattices. The method allows the 
systematic calculation of the main and corrective asymptotic terms to any order for large A^. 
These corrective terms are generally non-negligible as they (essentially) decay as powers of 
1/ In A^. However, we found that the main and first two corrective terms lead to reasonably 
good results when relatively small values of A^ are used (say, for A^ > 2''). In Sec. |T| we 
proposed a geometric meaning for the main and corrective terms: the main term would 
account for the volume of the set of visited sites if the exploration of the random walkers 
were compact, and the corrective terms just improve this rough estimate because, in the 
outer regions, the exploration performed by the (relatively few) random walkers that move 
there is really not compact, thus leading to the formation of a non-compact (a dendritic) 
external ring in the set of visited sites. We hope the above results and ideas could serve as 
a basis to gain insight into problems with interacting random walkers. 

The method developed here for calculating S'Ar(t) is also useful for evaluating other 
statistical quantities related to the diffusion of a set of independent random walkers. An 
example is the number SN+{t) of sites visited by A^ random walkers on an one-dimensional 
lattice along a given direction [|ll[]. It turns out that the moments (of arbitrary order) of 
Siy+(t) can be readily obtained through a slight modification of Eq. (pil). Another example 
is the first passage time ti,Ar(r) to a distance r of the first random walker of a set of A^. First 
passage times are relevant statistical quantities in the study of diffusion processes where 
the arrival of the first particles at a given site produces a significant effect (a "trigger" 
effect). These quantities have been calculated for one dimension |]3|,p^ (and for some classes 
of fractals |0]) but little is known for dimensions greater than one [Q]. The approximate 
compact form of the set of visited sites allows one to estimate the first passage time via 
the relation 5'Ar(ti,Ar(r)) ^ vor'^ [0, which means geometrically that we consider the region 
inside the hypersphere of radius r where a random walker has arrived by time ti^N{f) as 
completely visited (a compact exploration in the sense of de Gennes [|14|)- Results on 



SN+{t) and ti,Ar(r) obtained using the above ideas will be reported elsewhere. 

The function S'Ar(t) we have studied is indeed an important quantity concerning the 
diffusion of A^ independent random walkers but there are still many open questions in this 
problem. One can think, for example, of the absorption probability of the set of A^ random 
walkers on a lattice with a random distribution of point-like traps. This problem can be 
formulated in terms of the moments of the number of distinct sites visited by the set of A^ 
walkers. A prediction for the variance of the number of visited sites is a necessary requisite 
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to tackle this interesting problem as the first-order approximation based only on its first 
moment, i.e., on S]\f{t), seems to be very imprecise [§]. As no relationship is known for 
moments of order higher than one, the absorption problem remains unsolved. 

Finally, it should be pointed out that the expression for S'Ar(t) given in this paper can be 
extended to fractal media with some slight changes. We are currently running simulations 
for deterministic (Sierpinski gasket) and stochastic (percolation aggregate) fractals. Results 
for these substrates will be published elsewhere. 
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APPENDIX 



We will show in this Appendix how to get Eq. (|T^) from Eq. (0). Let us start by showing 
that the solution C,{z) of Eq. ([T6|) for z ^ has the form given in Eq. (pISl). For simplicity 
of notation we will write u = = 1 — exp(— 2;) and c = d/2. Hence, Eq. (|T^ takes the 
form 



h yulnu + InA + In ( 1 + ^ I = ln0 . (Al) 



n=l 



In the limit z — ^ 0, it is clear that m — *• and — > 0. This means that, as long as 
1/u ^ \ ln{A) I , the first term on the right-hand side of ( [Al| ) is the most divergent one so 
that, as a first approximation, we have 

u ^ --^ = X . (A2) 
ln0 

This first-order approximation was already obtained in Sec. f| [see below Eq. (0)]- ^ better 
approximation is achieved by writing u = x(l-|-e), with e an small quantity. The substitution 
of this expression in (^) yields 



2 ux X ux /ix 2 , hix'^ hix^e 
e — e H mx^ — mA^ e e H 1 h . . . = , (A3) 

c c c 2c c c 



where (|A2|) has been taken into account. This equation can be solved by writing e as 

e = enx" , (A4) 



n=l 



and inserting it in ([A 3]). We thus find the following values for en up to n = 2: 

= -iln(Aa;^), (A5) 

c 

e2 = \ In^ {Ax^) + 4 In {Ax^) - — . (A6) 
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Therefore 



-1/2 



X 



X 



(A7) 



n=0 



where 5q = 1 and 



5i 



In {Ax^) 
Yc ' 



(A8) 



(52 = ^ In^ (Ax'^) + 4/iln (Ax'') - 4c/ii 



The evaluation of the integral for Ki^{a) in (^) requires the expression of ^"'d^/dz as a 
function of From Eq. ( [A7|) and taking into account that d^/dz = {d^/dx){dx/dz) and 
dx/dz = [x'^/{c(f))]d(f)/dz, we find that 



e 



dz 



2c(f) dz 



X 



{l-a)/2 



1 + E ^" E k^r^^ In'" (^^") 



n=l m=l 



(A9) 



where the coefficients m = 0, . . . ,n for n = 1, 2 are 



(1) 



£.(2) 



c 

a — 1 

(a — 3)/ii ^ /i^ 



(2) 



2c c 
/i(2 — a) 

a{a - 4) + 3 



2 ' 



Let us use Ki^{a, z) to denote the integrand of Eq. (|I5D, i.e., 



Then, from Eq. ( |A9|) , 



K^ia^z) 



2c 



_g-Af^g-2^(l-Q)/2 



oo n 



1 + 5: X" ^ fc^") In- (A 



X'' 



n=l m=l 



(AlO) 



(All) 



(A12) 



Writing e"^ = 1 + 0{z), x = -(c/ln2)[l + 0{z/ In z)] and ln(y4x^) = InA - /xln(-lnz) + 
filnc + 0{z/ In z), Eq. ( |A12| ) becomes 



Kr,{a,z) = [l + 0{z)] 



-Nzi 



2c("+i)/2 



,ln^)("-i)/2 ^ ^ fc^)(-ln2)"ln'"(lnz) (A13) 



n=0 m=0 

where the coefficients fc^-* up to second order [n = 2) are 



11 



(1) 



(1) 



.(2) 



(2) 



(a - 1)- -/X, 
(1 -tt)2' 



/llC 



(3 - a)(l - a)— + /i(2 - + /i^ + ^ (a - 3) 
o 2 



(a-2)/x + (a-3)(l-a)- 



and a; = j + In A + filnc. Finally, inserting Eq. ( A13|) into Eq. ( All ) we get Eq. (19). It 
should be noted that we have approximated the factor l + 0{z) of Eq. ( [A13| ) by 1. This can 
be done safely because the contribution of the neglected terms to the asymptotic behavior of 
K^{a) decays as least as (In A^)'^"^-'^^/^/A^^, i.e., decays to zero faster than the contribution 
of the retained terms by (roughly) a factor [see Eqs. (p!9D-(p3D]. 
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TABLES 



TABLE L Parameters appearing in the asymptotic expression of S]\f{t), Eq. (^). The 
symbol dD refers to the d-dimensional simple hypercubic lattice. The parameter p is 
[2t(2L»7r)3/3]^''^p(0,l), where p(0,l) ~ L516386 §. 
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FIGURES 



FIG. 1. The integrand N^Vf-^ {dVt/di) of J^(l; 0, oo) versus ^ for the one-dimensional lattice 
and = 1, iV = 20 and N = 100. The solid lines correspond to the integrand when the exact value 
of Tt{S,) is used. The broken lines are obtained by using the first-order asymptotic approximation 
Vtii) ^ 1 - (2/7r)V2^-iexp[-^V2](l - ^"2). The filled circle [square] marks the value of for 
N = 20[N = 100] using p = 2 in Eq. (|). 

FIG. 2. A snapshot of the set of sites visited by = 1000 random walkers on the 
two-dimensional lattice. The visited sites are in white, the unvisited ones are in black and the 
internal gray points are the random walkers. The outer white circle is centered on the starting 
point of the random walkers and its radius is the maximum distance from that point reached by 
any walker at the time the snapshot was taken. The internal black circle is concentric with the 
former but its radius is the distance between the origin and the nearest unvisited site. 

FIG. 3. Four successive scaled snapshots of the set of sites visited by = 700 random walkers 
on the two-dimensional lattice for times (from left to right) t = 2000, t = 4000, t = 6000 and 
t = 8000. The second snapshot has been shrunk by the factor 1/^/2, the third by the factor l/\/3 
and the last by the factor 1/2. 

FIG. 4. S = [SN{t)/vo]'^/'^ / {4dDtln N) versus 1/lniV for, from top to bottom, dimension 1, 2 
and 3 and t = 200 (inside time regime II). We have used N = 2"^ with m = 3, • • • , 14 for d = 2,3, 
and m = 3, • • • , 30 for d = 1. The numerical results are plotted as filled circles and the broken 
[solid] lines correspond to the theoretical predictions for S^it) to first [second] order as given by 
Eq. (p4|). Notice that the approximation of order would be a horizontal line (not shown here) 
passing through l/d. The crosses correspond to the Sastry and Agmon result of Eq. ( |3^ ) with 
a = 1. The dotted lines correspond to the result of Larralde et al. given by Eq. (|l]) in which the 
corrected amplitude of the main term has been used (see p^). 
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